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ABSTRACT 

We investigate the variety in quasar UV spectra (AA1020-1600) with emphasis on the 
weak emission lines in the Lya forest region using principal component analysis (PCA). 

This paper is a continuation of Suzuki et al. (2005, Paper I), but with a different ap¬ 
proach. We use 50 smooth continuum fitted quasar spectra (0.14 < z < 1.04) taken by 
the Hubble Space Telescope ( HST) Faint Object Spectrograph. There are no broad ab¬ 
sorption line quasars included in these 50 spectra. The first, second and third principal 
component spectra (PCS) account for 63.4, 14.5 and 6.2% of the variance respectively, 
and the first seven PCS take 96.1% of the total variance. The first PCS carries Lya, 

Ly/3 and high ionization emission line features (O VI, N V, Si IV, C IV) that are sharp 
and strong. The second PCS has low ionization emission line features (Fe II, Fe III, 

Si II, C II) that are broad and rounded. Three weak emission lines in the Lya forest 
are identified as Fe II A1070.95, Fe II+Fe III A1123.17, and C IIP Al 175.88. Using the 
first two standardized PCS coefficients, we introduce five classifications: Class Zero and 
Classes I-IV. These classifications will guide us in finding the continuum level in the 
Lya forest. We show that the emission lines in the Lya forest become prominent for 
Classes III and IV. By actively using PCS, we can generate artificial quasar spectra that 
are useful in testing the detection of quasars, DLAs, and the continuum calibration. We 
provide 10,000 artificially generated spectra. We show that the power-law extrapolated 
continuum is inadequate to perform precise measurements of the mean flux in the Lya 
forest because of the weak emission lines and the extended tails of Lya and Ly/3/O VI 
emission lines. We introduce a correction factor 5F so that the true mean flux (F) can 
be related to {Fp ower _p a w) as measured using power-law continuum extrapolation by : 

(F) = (F Power _ Law )-5F. The correction factor SF ranges from 0.84 to 1.05 with a mean 
of 0.947 and a standard deviation of 0.031 for our 50 quasars. This result means that 
using power-law extrapolation we miss 5.3% of flux on average and we show that there 
are cases where we would miss 16% of flux. These corrections are essential in the study 
of the intergalactic medium at high redshift in order to achieve precise measurements 
of physical properties, cosmological parameters, and the reionization epoch. 
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1. Introduction 

It is important for the study of quasar absorption lines to understand the shape of the continua 
in the Lya forest from which we study the physical properties of the intergalactic medium (IGM) 
and extract cosmological parameters (Kirkman et al. 2003; Tytler et al. 2004a). It is the uncertainty 
of continuum fitting to the Lya forest that makes precise measurements difficult(Croft et al. 2002a; 
Suzuki et al. 2003; Jena et al. 2004). Thus, we wish to have a simple and objective quasar spectrum 
classification scheme that enables us to describe the global shape of the continuum as well as the 
individual emission line profiles. The principal component analysis (PCA), also known as Karhunen- 
Loeve expansion, is one of the best methods to carry out such classification. 

PCA enables one to summarize the information contained in a large data set, and it is widely 
used in many areas of astronomy (Cabanac et al. 2002, and references therein). Francis et al. 
(1992) applied PCA to quasar spectra using 232 quasar spectra (1.8 < z < 2.2 ; AA1150 — 2000) 
from the Large Bright Quasar Survey (LBQS; Hewett et al. 1995, 2001). They showed that the 
first three principal components account for 75% of the variance. Boroson & Green (1992) used 
87 low redshift quasar spectra ( z < 0.5) and showed the anticorrelation between the equivalent 
width of Fe II and [O III] and the correlation between the luminosity, the strength of He II A4686, 
and the slope index a ox . Boroson (2002) investigated the relation between the first two principal 
components and the physical properties such as black hole mass, luminosity, and radio activity. 
Shang et al. (2003) studied 22 low redshift UV and optical quasar spectra (z < 0.5 ; AA1171 —6607) 
and showed the relation between the first principal component and the Baldwin effect, which is the 
anticorrelation between the luminosity and the equivalent width of the C IV emission line. Yip 
et al. (2004) applied PCA to the 16,707 Sloan Digital Sky Survey quasar spectra (0.08 < z < 5.41 ; 
AA900 — 8000) and reported that the spectral classification depends on the redshift and luminosity, 
and that there is no compact set of eigenspectra that can describe the majority of variations. They 
also showed the relationship between eigencoefRcients and the Baldwin effect. 

In Suzuki et al. (2005, hereafter Paper I) we analyzed 50 continua fitted quasar spectra taken 
by the Hubble Space Telescope ( HST ) Faint Object Spectrograph. Since they are at low redshifts 
(0.14 < z < 1.04), and the Lycc forest lines are not so dense, we can see and correctly fit the 
continuum to the spectra. Using PCA we attempted to predict the continuum in the Lya forest 
where the continuum levels are hard to see because of the abundance of the IGM absorptions. 
Although we succeeded in predicting the shape of the weak emission lines in the Lya forest region, 
our prediction suffers systematic errors of 3 — 30%. This paper is a continuation of Paper I, but 
with a different approach. 

The goal of this paper is to explore the variety of quasar UV spectra in the following manner: 
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1. we clarify the PCA formulation in order to describe quasar UV spectra quantitatively (§3) using 
eigenspectra or the principal component spectra (PCS), 2. we introduce five classes of quasar UV 
spectra to help us understand the variety of quasar spectra qualitatively (§5), 3. we introduce the 
idea of artificial quasar spectra (§4) and 4. the mean flux correction factor 5F (§6), and 5. we 
report the identities of three weak emission lines in the Lya forest (Appendix). 


2. Data 

We use the same 50 HST FOS spectra from Paper I, and the detailed description is therein. 
Here we summarize the 50 HST spectra. These 50 quasar spectra are a subset of the 334 high 
resolution HST FOS spectra (R ~ 1300) collected and calibrated by Bechtold et al. (2002) which 
include all of the HST QSO Absorption Line Key Project’s data (Bahcall et al. 1993, 1996; Jannuzi 
et al. 1996). Bechtold et al. (2002) identified both intergalactic and interstellar medium lines and 
corrected for Galactic extinction. 

For each spectrum, we combined the individual exposures and remeasured the redshift using 
the peak of the Lyct emission line. Then we brought the spectrum to the rest frame and rebinned 
it into 0.5 A pixels. We masked the identified absorption lines and fitted a smooth Chebyshev 
polynomial curve. We normalized the flux by taking the average of 21 pixels around 1280 A where 
no emission line is seen. 

We chose the wavelength range from 1020 to 1600 A so that we can cover the Lya forest 
between the Ly j3 and the Lya emission lines while maximizing the number of spectra. We removed 
broad absorption line (BAL) quasars from the sample since we are interested in emission line profiles 
and the continuum shape. Lastly, we removed quasars with S/N < 10 per pixel because we cannot 
extract weak emission line features in low S/N spectra. The redshift range of the 50 spectra is from 
0.14 to 1.04 with a mean of 0.58. The average S/N is 19.5 per 0.5 A binned pixel. 


3. PCA 8z Principal Component Spectrum 
3.1. The PCA Formulation 

We express a quasar spectrum in Dirac’s “bra ket” form, \q%), which is commonly used in 
Quantum Mechanics and simplifies our description. We claim that any quasar spectrum, \qi) , is 
well represented by a reconstructed spectrum, |rj m ), which is a sum of the mean and the weighted 
m principal component spectra: 


\qt ) ~ I r i,m) ~ \fj) + y ' Cjj |£j) 
3 =1 


(1) 
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where i refers to a quasar, |//) is the mean quasar spectrum, | £j) is the jth principal component 
spectrum (PCS), and c^j is its weight. Unlike Quantum Mechanics, the weight, c ZJ , is not a complex 
number but a real number. We found the covariance and correlation matrix of the 50 quasar spectra 
in Paper I and by diagonalizing the covariance matrix, we can obtain the PCS. In practice, we 
found the PCS via a singular value decomposition (SVD) technique from the 50 quasar spectra. 
The recipes can be found in Francis et al. (1992) and Paper I. Since we designed PCS to be 
orthonormal: 


<&&> = Sij, 


where we define the inner product as: 


/■1600A 

(x\y) = / 0 x(A) y{ A) d\, 

J 1020A 


and we can obtain the weights Cj,- as follows: 


Cij = (Qi ~ lA £?)• 


( 2 ) 

( 3 ) 

( 4 ) 


We show the mean spectrum in Figure 1 and the first 10 PCS and the distribution of their 
weights in Figure 2 and 3. Although wavelength coverage, normalization, resolution, and the 
number of quasars are different, the general trend of PCS is similar to that of Francis et al. (1992) 
with the exception of the third PCS. The reason for the exception is that the third PCS of Francis 
et al. (1992) takes into account BAL features, while we don’t have such features since we removed 
BAL quasars. Our third PCS has sharp emission line features for Ly/3/O VI, the core of Lyct, and 
C IV emission lines, but broad negative features around the Lyct emission line. The continua of 
the third, fourth and fifth PCS have a slope, and the fifth PCS includes the asymmetric feature for 
the strong emission lines. The fifth PCS is similar to the second PCS in that P-Cygni profiles are 
seen but they are very broad, and there are no low ionization emission line features blueward of 
the Lya emission line. The sixth and seventh PCS carry some information on low ionization weak 
emission lines, but the spectrum features are getting noisy. The eighth and higher order PCS have 
high frequency wiggles which no longer correspond to any physical emission lines and are probably 
due to fitting errors. We used continua fitted spectra which are free from photon noise, but still 
they are likely to suffer fitting errors of at least a few percent as we can see in Figures 6-10. The 
mean spectrum and the first 10 PCS are available on-line from the Paper I. Next, let us discuss the 
contribution from each PCS quantitatively. 


3.2. Quantitative Assessment of PCA Reconstruction 

We make use of the residual variance to assess the goodness of PCA reconstruction. The 
residual variance of a reconstructed quasar spectrum with m components is a sum of the squares 
of the difference between the quasar spectrum, | q^), and the reconstructed spectrum, | ri <m ): 
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N 

(Qi ~ r ijTn \qi - r iiTn ) = ^ c?-, (5) 

i=m -\-1 

where N is the total number of independent PCS; either the number of quasar spectra or the total 
number of pixels. In our case, N is the number of quasar spectra (50) since it is smaller than the 
total number of pixels (1167). For the overall contribution from mth PCS, we define the residual 
variance fraction, i(j ), as follows: 


fti) = 


N 

£■ 

i= 1 


N 

£■ 

** 7=1 


( 6 ) 




We redefine the square root of the eigenvalue of the jth PCS as follows: 


A, 


1 

N- 1 


N 



( 7 ) 


Then we wish to use A j to describe the probability distribution function (PDF) of the PCS 
coefficients. The PDF of the weights is not necessarily a Gaussian, but as shown in Figures 2 and 
3, the PDF can be well represented by a Gaussian. We also note that by design: 


N 

'y ^ c i j = 0 (s) 

i =1 

for any j. Thus, the PDF of weights is characterized by just one parameter, A j. The probability of 
having a weight, c t j, in an interval, — xq < Cij < xq can be expressed as: 


P(-xq < Cij < 





( 9 ) 


Naturally, it would be convenient if we standardize the weight, Cjj, by A j. We can then rewrite 
a quasar spectrum as: 

m 

I Qi) ~ I r i,m) = ll) y ^ AjtTjj |^j) (10) 

3 = 1 

where the cjjj are the standardized PCS coefficients which represent the deviation from the mean 
spectrum of the jth PCS of quasar i. The PDF of the a tJ is a normal distribution, so we can 
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immediately tell how far and how different the quasar spectrum would be from the mean spectrum. 
This standardization of the weights simplifies the discussion of the variety of quasar spectra in §5. 

Another advantage of using A j is that we can simplify the residual variance fraction f(j): 

m = ~n~' f 11 ) 

3 =1 

The values of f(j) are listed in Table 1. The first, second and third PCS account for 63.4%, 
14.5% and 6.2% of the variance respectively. In total, the first three PCS account for 84.3% of 
the variance of the 50 quasars in our sample. These fractions depend on the normalization and 
the wavelength coverage. In the literature, Francis et al. (1992) report that the first three PCS 
account for 75% (AA1150 — 2000), and Shang et al. (2003) show that the first three PCS carry 
80% (AA1171 — 2100) of the variance. Both groups normalized flux by the mean flux, while this 
paper normalizes by a flux value near 1280 A. Our value of 84.3% is slightly higher than the above 
numbers probably because we removed the BAL quasars, which are certainly a source of variance. 
In addition, we used fitted smooth continua to the Lya forest, while they used the observed raw 
Lya forest which obviously contains a large variance (Tytler et al. 2004a). As shown in Table 1, the 
contribution from each PCS component to the variance rapidly decreases with order m. It becomes 
less than 1% after the 8th PCS and then stays the same. With seven PCS components, 96% of the 
variance has already been accounted for. As is seen in the PCS features in Figure 3, the remaining 
4% of the variance is probably due to fitting errors. 


4. Artificial Spectra 

We can use PCS to generate artificial spectra. Artificial spectra can be useful in testing 
the detection of quasars and DLAs, in flux calibration, in continuum fitting, and in cosmological 
simulations. By assigning PCS coefficients randomly from known PDFs, we can generate artificial 
quasar spectra. As we have discussed in §3, the PDF of the jth PCS coefficient is well represented 
by a Gaussian with a mean of 0, and a standard deviation of A j. If we then sum up the PCS with 
these coefficients (equation 10), we can create a set of artificial quasar spectra. 

Noiseless quasar spectra are of great use in IGM studies, since at high redshift, it is difficult to 
see the continuum in the Lya forest. Even at redshift 2, pixels in the Lya forest hardly reach the 
continuum level with FWHM=250km/s (Tytler et al. 2004b). Artificial quasar spectra can thus be 
useful to predict the shape of continuum in the Lya forest (Paper I) and to calibrate continuum 
fitting to the Lya forest (Tytler et al. 2004a,b). They would be also useful to test the detection 
limit of the high redshift quasar survey since the Lya emission can possibly boost the brightness by 
0.15 magnitude. We have generated 10,000 artificial quasar spectra using the first seven PCS for 
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a more realistic representation of the quasar spectra. We concluded that the features seen in PCS 
greater than eighth are noise, and we did not include higher order PCS. We will provide artificially 
generated spectra to the community upon request. 


5. PCA Classification 
5.1. Introduction of Five Classes 

This paper attempts to classify quasar spectra quantitatively using our standardized PCS 
coefficients, aij. We note that there is no discrete classification of quasar spectra and that they 
vary continuously. However, this classification will help us in fitting the continua to the Lya forest 
spectrum. We use the first two PCS coefficients to demonstrate the variety of emission lines and 
continua. As we have seen in §3 and Table 1, the first two PCS coefficients account for 77.9% of the 
variance and represent the overall shape of the quasar spectrum. We introduce polar coordinates 
as follows: 


r i = \j ( 12 ) 

tan Oi = —— (13) 

where r* represents the deviation from the mean spectrum, and the angle 6{ tells us about the 
profiles of the emission lines. We divide the an vs. a^ diagram into five zones and introduce five 
classes of spectral types that enable us to discuss the shape of continua qualitatively. The main 
goal of this classification is to differentiate the families of quasar spectra. First, we define Class 
Zero for those close to the mean spectrum in shape. We define Class I-IV corresponding to the 
quadrant I-IV in the an vs. aa diagram, shown in Figure 4. The probability of r < ro is: 

P(r < r 0 ) = 


Now, we wish to design the fraction of the five classes to be equal, namely 0.2 each. We find 
ro = 0.668 gives P(r < ro) = 0.2, so we define Class Zero for a quasar spectrum that has r < 
0.668. For quasar spectra with r > 0.668, we define the classes / through IV corresponding to the 
quadrants first through fourth on the a t \ vs. a^ diagram (Figure 4). 

Our standardized PCS coefficients are plotted on the an vs. a *2 diagram in Figure 4 where 
the dotted circle is the circle with radius ro = 0.668. The small numbers noted beside the points 
are the quasar identification numbers i that are listed in Table 2. 


/ 

Jo 


r 0 _ r 2 
re 2 dr 


1 — e 


(14) 

(15) 



5.2. Demonstration of Five Classes 


We show the mean spectrum in Figure 1. By definition the mean spectrum \/i) has r=0, and 
naturally it belongs to Class Zero. In Figure 5 we show the artificially generated four classes, I-IV , 
of the quasar spectra to illustrate their typical spectral shape. They are the sum of the mean and 
the first two PCS with an = ±1 and a^ = ±1. The generated four spectra of Class I-IV have 
angles 9: 45°, 135°, 225° and 315° respectively, and they all have r = \/2. The four spectra are 
plotted on the same scale in Figure 5 so that we can see the contrast of the emission lines with the 
continuum in a uniform manner. 

In Figure 6-10, we show three observed spectra from each class where we intentionally chose the 
extreme cases for Class I-IV. Quasars are plotted at rest frame wavelengths with the luminosities 
that are calculated by using the cosmological parameters from the first year WMAP observation 
(h=0.71, fl m =0.27, 14 a= 0.73 ; Spergel et al. 2003), and a flat universe is assumed. The smoothed 
line on the spectrum is the fitted continuum, and the solid straight line is the power-law continuum 
fit. The vertical dotted lines show the wavelengths of emission lines in the Lya forest and the low 
ionization lines redward of the Lyct emission line. The quasar numbering in Figure 6-10 is the same 
as in Figure 4 so that we can visualize where the quasar spectrum is in the an vs. < 7*2 diagram. In 
Table 2, quasars are sorted by classes, and the equivalent width of the emission lines are listed. 


5.3. The Characteristics of the Five Classes 

The characteristics of the first two PCS directly reflect on the five classifications, thus let us 
take a close look at the first two PCS in Figure 2. The first PCS carries the sharp and strong lines 
- Lya, Ly/3 and high ionization emission line features (O VI, N V, Si IV, C IV). The second PCS 
has low ionization emission line features: Fe II and Fe III in blueward of the Lya emission line, 
Si II and C II redward. Their profiles are broad and rounded. In the second PCS, the flux values of 
low ionization emission lines and the strong Lya, C IV emission lines have opposite sign, meaning 
they are anticorrelated. In addition to that, Lya and C IV emission lines have P-Cygni profiles 
which introduces asymmetric profiles to these emission lines. 

By definition, the first two PCS engage the correlation between emission lines. Since we are 
particularly interested in the profiles of emission lines in the Lya forest, let us look at low ionization 
lines first. If a quasar shows prominent low ionization lines redward of the Lya emission line (Si II 
AA1260,1304, C II A1334), it should have a negative second PCS coefficient, and we should expect 
to have prominent Fe II A1070 and Fe III A1123 in the Lya forest. Thus such a quasar should 
belong to either Class III or IV. As a consequence, these two classes have the largest equivalent 
widths of these low ionization lines among the five classes as seen in Table 2. 

If another quasar has sharp and strong Lya, Ly/3 and high ionization lines (N V, Si IV, C IV), 
it should have a positive first PCS coefficient and belong to Class I or IV. The normalized flux 
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of the Lya emission peak and the ratio of the Lya and N V peak flux are the highest for Class I 
among the five classes. We can differentiate Class I and II, or Class III and IV by combining 
the above characteristics. As we expect, the diagonal classes have the opposite characteristics. For 
example, Class I has sharp and high ionization lines while Class III has broad and rounded low 
emission lines. 

In practice, the key point of finding the continuum level in the Lya forest is to seek the low 
ionization lines (Si II AA1260,1304, C II A1334) and their profiles redward of the Lya emission. If 
we see them, we should expect to have similar profiles of Fe II A1070 and Fe II A1123 lines in the 
Lya forest. If we do not see them, we can expect the continuum to be flat in the Lya forest, and 
the power-law extrapolation from the redward of the Lya emission to be a good approximation. 
We will discuss the accuracy of the power-law extrapolation in the next section. 


6. Mean Flux (F) and Flux Decrement Da 

In performing a precise measurement of the IGM flux decrement D a, it is essential to include 
the effect of weak emission lines found in the Appendix A. This contribution is not negligible, and 
will be quantitatively assessed. 


6.1. A Brief History of the Flux Decrement D^ 

Gunn & Peterson (1965) predicted a flux decrement due to foreground neutral hydrogen in 
the IGM, namely the Lya forest. Oke & Korycansky (1982) first defined and measured the flux 
decrement, Da- Schneider et al. (1991) introduced the Lya forest wavelength interval for Da as 
AA1050 — 1170, and it has been widely used since (Zuo & Lu 1993; Kennefick et al. 1995; Spinrad 
et al. 1998). Madau (1995) and McDonald & Miralda-Escude (2001) have used D a to estimate the 
UV background. The flux decrement of the high redshift IGM probes the reionization epoch of the 
universe (Loeb &; Barkana 2001). The current estimate of the reionization epoch from the IGM is 
around £ ~ 6 — 7 (Becker et al. 2001; Djorgovski et al. 2001; Fan et al. 2003), while the first year 
WMAP satellite estimates z ~ 20 (Spergel et al. 2003). The discrepancy is yet to be resolved or 
explained (Cen 2003). 

A precise measurement of the flux decrement, D^, is of great importance for studies of the 
IGM (Rauch 1998) because it is very sensitive to the cosmological parameters, erg (the amplitude of 
the mass power spectrum) and Q,\, as well as to the UV background intensity (Tytler et al. 2004b; 
Jena et al. 2004). However, it is this sensitivity that makes the measurement a major source 
of error (Hui et al. 1999; Croft et al. 2002b). 
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6.2. The Mean Flux Correction Factor SF 

The flux decrement is defined as: 


Da - 1 - I " 70A 


/ 

JlQ 


I 1050 A /a(A; Continuum) 
Thus, what we are measuring is the mean flux: 


(F)= I 

J li 


1170A 


/a(A; Observed) 
1050 A /a(A; Continuum) 


dX 


(16) 


(17) 


However, the unabsorbed continuum level is not seen in the Lya forest and the power-law 
extrapolation from redward of Lya emission has been used as a continuum. In fact, what is 
reported in the literature as the mean flux is: 


(• F P , 


ower—Law) — 


L 


1170 A 


/a(A; Observed) 


-dX 


(18) 


I 1050 A /a (A; Power-Law) 
which is not exactly the same as (F) since the power-law is a crude approximation of the continuum 
in the Lya forest as we have seen in §5, Appendix A, and Figures 6-10. We wish to introduce a 
correction factor 5F : 

rihOA Power—Law) 

' 1050A 


SF = 


-dX 


/a(A; Continuum) 

so that we can estimate the true mean flux ( F ) from the reported mean flux ( Fp ower _L aw )'■ 


(19) 


(F) = (F Power _ Law ) • SF. (20) 

To calculate SF, we need to find the power-law extrapolation from the redward Lya emission. 
Since our wavelength range is limited and not as large as other survey data, it is not as easy to 
extrapolate. Moreover, we have a series of emission lines, and it is hard to define a continuum 
level with no emission lines. For example, as shown in Figure 9, Class III quasars have emission 
lines throughout this wavelength range. However, we have a fitted continuum in the Lya forest, 
and we take advantage of it. We choose two points and find a power-law fit which runs through 
these two points. We choose one from blueward (A1100) and the other from redward (A1450) of 
the Lya emission. Then, we can have an inter- and extra-polated power-law continuum and its 
exponent a v , for f v oc v av . The average of a u is -0.854 with the standard deviation of 0.507. The 
power-law continua show that they are all sensible first order approximations which well represent 
the continua in the Lya forest. The power-law continua are shown in Figure 6-10. 
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The calculated SF is listed in Table 2 and the average of 5F is 0.947, with a standard deviation 
of 0.031. This result means that the power-law approximation misses 5.3% of flux from the quasar in 
the wavelength range, and proves that the power-law approximation is inadequate to perform a 
D a measurement that attempts 1% accuracy (Tytler et al. 2004a; Jena et al. 2004). The distribution 
of SF is shown in Figure 11, and it is not a Gaussian. The PDF of 5F is asymmetric and has long 
tail toward small SF value. 

There are two major reasons for the missing flux. The first reason is that the Dyi wavelength 
range is still in the tail of the prominent Ly/3/O VI and Lya emission lines. For example, the 
blueward tail of the Lya emission starts near A1160 which is 10A below the A1170 upper limit of 
Da wavelength range. The effect from the tails of Ly (3 and Lya emission is common for all of the 
quasars as we can see in Figure 6-10. This contribution is about 4-5%, and we always miss this 
fraction of the flux, meaning 5F is always less than unity. 

The second reason is the contribution from the weak emission lines in the Lya forest. The 
continuum level at the emission lines is naturally above the power-law extrapolation, therefore, we 
would always expect to miss flux from the emission lines, making 5F always less than unity. The 
low ionization emission lines in the Lya forest are prominent for Class III & Class IV quasars. Four 
quasars in Class III Sz Class IV have 5F less than 0.9, meaning we miss more than 10% of the flux 
if we use power-law extrapolation. The quasar which has the smallest 5F, 0.86, is Q1544+4855. 
This quasar is shown on the top panel of Figure 10 and the power-law continuum fit looks sensible. 
Together with the tails of Ly/1 and Lya emissions, the low ionization emission lines in the Lya 
forest (Fe II, Fe III) contribute 14% of flux which is significant and should not be neglected. 

We expect that we need to apply at least a SF = 0.947 correction to the D a measurements 
in the literature. The discrepancy between the past Da measurements and Bernardi et al. (2003) 
using the weak emission profile fitting method is shown in Tytler et al. (2004a, Figure 22). The 
disagreement is approximately 5% at redshift 2 < z < 3, and in terms of the mean flux, the power- 
law fitted values (Press et al. 1993; Steidel Sz Sargent 1987) are always above that of profile fitted 
ones (Bernardi et al. 2003). The correction factor, SF = 0.947, explains this disagreement well. 

However, we expect that SF changes with redshift, and it is crucial to include the effect from 
the weak emission lines to investigate the reionization epoch. Known as the Baldwin effect (Baldwin 
1977), the emission profile of lines, such as C IV, and the luminosity of the quasar are correlated. 
Because of the anticorrelation between equivalent width of C IV and the luminosity, we expect that 
Class III quasars to be the brightest since they have the smallest C IV equivalent width. Due to the 
selection effect, we would expect to observe the brightest quasars at high redshifts. Therefore, the 
fraction of classes for the observed quasars should change with redshift. In fact, the highest redshift 
quasars reported by Fan et al. (2001), Becker et al. (2001) and Djorgovski et al. (2001) probably 
belong to Class III , because they all show weak Si II A1304 emission line implying they have weak 
emission lines in the Lya forest regain. J104433.04-012502.2 (z=5.80) and J08643.85+005453.3 
(z=5.82) definitely belong to Class III since they have broad Lya, N V, and Si II emission lines 
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and low Lya/N V emission intensity ratio. The mean flux correction 5F of Class III quasars are 
in the range of 0.91-0.97, meaning the power-law extrapolation misses 3-9% of the flux for Class 
III quasars. We note that the reported 1 -a error of the residual mean flux at redshift z=5.75 by 
Becker et al. (2001) is 0.03. The contribution from 5F is bigger than their claimed error, and it is 
systematic. This correction would bring the observed mean flux down by 3-9% percent and would 
bring up by the same fraction. Therefore, it is crucial to take into account the 5F correction 
in order to investigate the reionization epoch. 


7. Summary 

We analyzed the wide variety of the emission line profiles in the Lya forest both in a quanti¬ 
tative and qualitative way. We used PC A to describe the variety of quasar spectra, and we found 
that 1161 pixels of data (AA1020 — 1600 with 0.5 A binning) can be summarized by the primary 
seven PCS coefficients because the pixels are not independent but are strongly correlated with each 
other. We presented, for the first time, the idea of generating artificial quasar spectra. Our artificial 
quasar spectra should be useful in testing the detections, in calibrations, and in simulations. We 
introduced five classes to differentiate the families of quasar spectra, and showed how the classi¬ 
fication can guide us to find the continuum level in the Lya forest. It is essential to account for 
emission line features in the Lya forest to perform a precise measurement of the mean flux in order 
to probe cosmological parameters, the UV background, and the reionization epoch, otherwise, on 
average, the commonly used power-law extrapolation continuum misses 5.3% of the flux, and we 
have cases when it misses up to 14% of the flux. 

To investigate the high redshift Lya forest, we showed the need to account for the emission 
lines in the Lya forest. An emission line profile continuum fitting method by Bernardi et al. (2003), 
or improvement of the PC A method in Paper I would be useful for a large data set such as the 
Sloan Digital Sky Survey. If we can study the redshift evolution of the quasar spectra and if we can 
estimate the constituents of classes at a certain redshift, we would be able to estimate the mean 
flux statistically using the mean flux correction factor 5F. For precision cosmology, the formalized 
we presented here should play an important role. 
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STScI and by NSF grant AST-0098731. 


-13- 


A. Appendix : Weak Emission Lines in the Lya Forest 

It has been suggested that there exist weak emission lines in the Lya forest (Zheng et al. 1997; 
Telfer et al. 2002; Bernardi et al. 2003; Scott et al. 2004), however, their identities, strengths, and 
line profiles are not well understood. More importantly, we wish to know how they are correlated 
with other emission lines so that we can predict the strength and profiles of the emission lines 
in the Lya forest. We attempt to find the identities of the three weak emission lines in the 
Lya forest reported by Tytler et al. (2004b). In this paper, we use their measured wavelengths : 
A = 1070.95, A = 1123.17, A = 1175.88 A. 

A.l. Line at A = 1070.95 A : Fe II 

Zheng et al. (1997) and Vanden Berk et al. (2001) identified this line as Ar I A1066.66, but the 
contribution from Ar I cannot be this large. The Ar I line has another transition at A1048.22 whose 
transition probability (H ? ;fc = 4.94x lO 8 ^” 1 ) is stronger than that of A1066.66 (H 7 ;fe = 1.30 x lO 8 ^^ 1 ). 
However, there is no clear sign of the A1048.22 line feature in 50 HST spectra or 79 quasar spectra 
in Tytler et al. (2004a). In addition, the A1070 line has a broad and asymmetric profile which 
suggests that this line is a blend of multiple lines. Thus, A1070 is not likely to be a single Ar I line. 

Telfer et al. (2002) labeled this A1070 as N II + He II line, and mentioned S IV as a possible 
candidate. Scott et al. (2004) identify four lines in the proximity of this wavelength : S IV doublet 
(AA1062, 1073) and N II+He II+Ar I (A1084). As we have seen in §5, the A1070 li n e correlates 
with low ionization lines such as Si II and C II. S IV does not fit into this category. N II seems to 
be a good candidate, but the N II lines peaks around 1085 A which is 15 A away from what we see. 
It is unlikely that we have 15 A of wavelength error. For the same reason, He II (1085A), a high 
ionization line, is not likely to be the dominant contributor. However, it is reasonable to expect 
to have an He II line since He II A1640 is often seen in quasar spectra. This A1070 line is seen in 
the Q1009+2956 and Q1243+3047 spectra for which we have high S/N Keck HIRES spectra with 
a FWHM=0.0285 A resolution at this wavelength (Buries & Tytler 1998; Kirkman et al. 2003). If 
A1070 is comprised of the lines suggested by Scott et al. (2004), we would be able to resolve the 
individual lines. However, none of the individual emission line are resolved. This fact implies that 
this emission line is comprised of numerous weak lines and that they are probably low ionization 
lines because of the good correlation with other low ionization lines. 

Fe II suits such a description, and in fact, Fe II has a series of UV transitions around this 
wavelength range: AA1060-1080. However, we don’t see other expected Fe II emission lines. If this 
A1070 is Fe II, we would expect to see the Fe II UV10 multiplet around A1144 , which is supposed 
to be stronger than the A1070 line. But no sign of an emission line is seen at that wavelength. 
Therefore, Fe II identification may be wrong or there may be a mechanism that we are not aware 
of that prevents the expected A1144 emission line. 
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A.2. A = 1123.17 A : Fe II + Fe III 

Telfer et al. (2002); Vanden Berk et al. (2001) identify A1123 as Fe III. The f-value weighted 
Fe III UV1 multiplet has a wavelength at 1126.39 A which is very close to what is observed by 
Tytler et al. (2004a). The distribution of the multiplet lines matches the broad A1123 feature found 
in quasars. There are no other major resonance lines in this wavelength range except C I. In the 
wavelength range AA1114 — 1200, C I has a series of lines, and there exists stronger lines redward 
of the Lya emission line: AA1115 — 1193, AA1277— 1280. However, we do not see these redward C I 
emission lines, and there is no sign of correlation between A1123 and these possible lines (Paper I). 
Thus, A1123 is probably Fe III. In addition to Fe III, Fe II also has UV11-14 multiplets around this 
wavelength: AA1121 — 1133. Given the fact that this A1123 line is well correlated with A1070 line, 
it is reasonable to expect to see Fe II lines here as well, if the identification of A1070 is Fe II. 


A.3. A = 1175.88 A : CIII* 

Telfer et al. (2002); Vanden Berk et al. (2001) identified A1175 as C IIP, although as shown 
in Table 2, we do not have a clear detection of this line because of it being a weak and narrow 
feature. We can see the A1175 line in the HST composite spectrum (Telfer et al. 2002), the SDSS 
composite spectrum (Vanden Berk et al. 2001), and in 11 out of 79 quasar spectra in Tytler et al. 
(2004a). Therefore, we are confident that this A1175 is a real emission line. The /—value weighted 
wavelength of the C IIP line is Al 175.5289 which matches well the observed wavelength. There is 
no major resonance line at this wavelength. 


A.4. Other Possible Emission Lines in the Lya Forest 

Tytler et al. (2004b) and Telfer et al. (2002) reported observing the Si II A1195 line in their 
spectra. We do not have any clear detection of the Si II A1195 line in the 50 quasar spectra. Since 
other Si II lines, AA1265, 1304, are clearly seen redward of the Lya emission line, it is plausible to 
expect the Si II A1195 line. 

Si III has a transition at A1206.50. Since we see both Si II and Si IV redward of Lya emission, 
it is natural to expect to see Si III. However, no detection is reported in the literature, and we do 
not have any positive detection in the 50 quasar spectra. The Si III A1206 emission line is probably 
too weak or possibly overwhelmed by broad Lya emission, only 10 A away, whose equivalent width 
is often greater than lOOA. 
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Rest Wavelength (A) 

Fig. 1.— The mean spectrum of 50 HST quasar spectra: The spectrum is normalized near 1280 
A. The wavelengths are taken from Morton (1991) except Fe II, Fe III and CIII* lines that are 
observed wavelenghs from Tytler et al. (2004b). The tick marks shown below the spectrum are 
the wavelengths of the Fe II and Fe III multiplet. The dotted line is the power-law continuum 
apporximation. Note the emission lines do exit in the Lycc wavelength region, We also note the 
wavelength separation of Si IV doublet at A1400 is relatively large, and makes the line profile broad. 
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Fig. 2.— The first five PCS 1-5 are shown in the left panels with the distribution of PCS coefficients 
Cij on the right. The first PCS takes the Lya, Ly/3 and the high ionization line features (O VI, 
Si IV, C IV) while the second PCS account for the low ionization lines (Fe II, Fe III, Si II, C II). We 
expect to have the low ionization lines in the Lya forest when a quasar has the negative second PCS 
coefficient. We note that PCS are designed to be orthonormal. The quasar spectra are normalized 
near 1280 A and every PCS has zero value near 1280 A. 
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Fig. 3.— The same as Figure 2 but the second five PCS 6-10. From 8th and higher PCS, the 
features in PCS become noisy and no longer have physical correspondence. They are probably 
fitting errors. 
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Fig. 4.— The distribution of standardized first two PCS coefficients for 50 quasars. The small 
number besides each point is quasar number i which is listed in Table 2. For quasars with open 
circles, we have quasar spectra shown in Figure 6-10. We intentionally chose extreme cases to show 
a wide variety of quasar spectra. The the dotted circle has the radius of 0.668. We divide this 
plane into the five zones and define the five classes. We define Class Zero as the quasars which 
have radius less than or equal to 0.668 and Class I-IV as the ones which have radius greater than 
0.668 and are in quadrants I-IV respectively. 
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Fig. 5.— The illustration of four classes. We artificially generated these four spectra using the first 
two PCS. We chose and ±< 7*2 to generate the four spectra, and equation is given in the legend 
where |/r) is the mean spectrum, Ai and A 2 are the square root of the first two eigenvalues, and 
|£i) and |^} are the first two PCS. The wavelengths of low ionization lines are shown as vertical 
dotted lines. Note that the four spectra are plotted in the same scale, and the emission line peak 
contrast with the continuum characterizes the classification. Class I and II do not have any weak 
emission lines and the power-law continuum fit (smooth dotted line) is a good approximation of 
the continuum in the Lya forest while Class III and IV show prominent low emission lines. 
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Fig. 6.— Class Zero ( r < 0.668) : The smooth line on the spectrum is the fitted continuum, and 
the straight solid line is the power-law continuum fit. The line at the bottom is 1 -a error of the 
spectrum. The vertical dotted lines show the wavelengths of emission lines in the Lya forest and 
low ionization lines redward of the Lya emission line. The luminosity, L^, is measured at A1280 
where we normalized the spectra. The emission lines in the Lya forest are barely seen in the three 
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Fig. 7.— Class I ( r > 0.668, 0 < 0 < 90°): The Lya, Ly /3 and high ionization emission lines are all 
very sharp and high contrast with the continuum. The low ionization lines, Si II at AA1260,1304, are 
barely seen. We see no clear sign of emission lines in the Lya forest, and the power-law continuum 
is a good approximation. The contrast of Lya emission peak with the continuum is the highest 
among the five classes. In Q2344+0914, there is a huge spike at 1030 A. It is due to an emission 
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Fig. 8.— Class II ( r > 0.668, 90° < 9 < 180°): The Ly a emission peak has a moderate contrast, 
4-6, with the continuum, and the emission peak ratio of Lya and N V is 2-3. The Lya emission 
has long tails, and the blueward tail of the fitted continuum does not meet with the power-law 
continuum until 1120 A, which is 50 A below the wavelength definition, 1170 A. In an extreme 
case, Q0349-1438 (r = 2.345), the tail of Lya emission wings last about 100 A. Ly/3/O VI, Lya, 
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Fig. 9.— Class III (r > 0.668, 180° < 9 < 270° ): The emission line profiles are all broad and 
rounded. Fe II and Fe III lines are clearly seen in the Lya forest. The contrast of the Lya emission 
line peak with the continuum is the lowest, 2-4, among the five classes. The ratio of Lya emission 
peak to N V is also the lowest: 1-2. The C IV profile is asymmetric probably because of the 
contribution from the Si II A1526 hidden inside the blueward tail of the C IV line. 
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Fig. 10.— Class IV (r > 0.668, 270° < 6 < 360°): The high contrast with the continuum and 
sharp emission profile are the characteristics of this class. Q1544+4855 is an extreme case that 
has r=3.099. Since it has 8 = 278.5, it is very close to Class III, but it has the characteristics of 
Class IV: high Lya emission line peak contrast and the very sharp line profiles. Unlike a Class III 
quasar, the peaks of low ionization line profile are very sharp and not rounded. Q1322+6557 has a 
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Fig. 11.- The distribution of the mean flux correction factor 5F. The mean is 0.947 with the 
standard deviation of 0.031. The distribution is asymmetric with long tail toward small 5F values. 
Class III and IV quasars have prominent emission lines in the Lycc forest, and tend to have small 
5F values. Two extreme spectra for low 5F are shown in Figure 9 for Q0959+6927 and Figure 10 
for Q1544+4855. 
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Table 1. Eigenvalue and Residual Variance Fraction 


Component : j 

i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Eigenvalue 

57.204 

12.985 

5.528 

4.614 

2.516 

2.189 

1.293 

0.605 

0.541 

0.453 

V 

7.563 

3.604 

2.351 

2.148 

1.586 

1.479 

1.137 

0.778 

0.735 

0.673 

Residual Variance Fraction f (j) 

0.637 

0.145 

0.062 

0.051 

0.028 

0.024 

0.014 

0.007 

0.006 

0.005 

Cummulative Residual Variance Fraction 

0.637 

0.781 

0.843 

0.894 

0.922 

0.947 

0.961 

0.968 

0.974 

0.979 


Note. 


Cummulative residual variance fracion is a simple sum of residual variance fraction upto j th PCS. 






Table 2 


The mean flux correction factor 5F and equivalent width of emission lines 


i 

Quasar 

< 5F 

Ly / 3/OVI 

1025 

(A) 

Fe II 

1071 

(A) 

Fe III 

1123 

(A) 

C III * 

1176 

(A) 

Lyo : 

1216 

(A) 

N V 

1240 

(A) 

Si II 
1263 

(A) 

Si II 

1307 

(A) 

C II 

1335 

(A) 

Si IV 

1397 

(A) 

C IV 

1549 

(A) 

Class Zero 

13 

Q0947+3940 

0.972 

8.8 

1.4 

0.0 

0.0 

102.9 

0.7 

0.5 

0.0 

0.0 

13.2 

59.8 

25 

Q1229-0207 

0.953 

9.5 

0.9 

0.2 

0.0 

87.7 

0.4 

0.1 

1.8 

0.0 

6.9 

47.4 

34 

Q1424-1150 

0.903 

9.5 

1.1 

0.0 

0.0 

131.9 

0.6 

0.0 

0.0 

0.0 

8.1 

50.4 

37 

Q1538+4745 

0.926 

8.2 

0.4 

0.9 

0.0 

132.9 

0.0 

0.0 

0.4 

0.0 

6.8 

43.4 

48 

Q2340-0339 

0.963 

5.1 

0.3 

0.0 

0.0 

83.5 

0.7 

0.0 

0.0 

0.0 

5.7 

41.5 


Average 

0.944 

8.2 

0.8 

0.2 

0.0 

107.8 

0.5 

0.1 

0.4 

0.0 

8.1 

48.5 


STD 

0.028 

1.8 

0.5 

0.4 

0.0 

23.6 

0.3 

0.2 

0.8 

0.0 

2.9 

7.2 

Class I 

1 

Q0003+1553 

0.940 

8.6 

0.1 

0.0 

0.0 

115.7 

0.0 

0.0 

2.2 

0.0 

6.5 

55.2 

8 

Q0439-4319 

0.970 

13.3 

0.0 

0.0 

0.0 

119.2 

1.1 

0.0 

1.2 

0.0 

7.3 

58.0 

11 

Q0637-7513 

0.980 

15.6 

2.5 

0.0 

0.0 

87.2 

2.2 

0.3 

0.0 

0.0 

8.6 

38.9 

14 

Q0953+4129 

0.937 

12.1 

1.9 

0.3 

0.0 

124.2 

1.4 

0.3 

0.0 

0.0 

10.6 

63.6 

18 

Q1007+4147 

0.943 

11.6 

2.0 

0.1 

0.1 

145.4 

1.9 

0.0 

2.0 

0.3 

10.4 

65.5 

22 

Q1137+6604 

0.939 

10.8 

1.8 

0.2 

0.0 

99.2 

0.6 

0.2 

0.6 

0.0 

7.3 

45.1 

24 

Q1216+0655 

0.959 

9.9 

0.9 

0.0 

0.0 

124.7 

0.0 

0.0 

0.6 

0.0 

8.8 

52.4 

32 

Q1354+1933 

0.974 

10.3 

0.2 

0.1 

0.1 

134.6 

0.5 

0.5 

2.0 

1.1 

7.9 

65.5 

39 

Q1622+2352 

0.971 

9.3 

0.7 

0.0 

0.1 

111.3 

0.2 

0.0 

0.0 

0.0 

7.7 

65.2 

41 

Q1821+6419 

0.986 

13.6 

0.3 

0.0 

0.0 

115.3 

0.2 

0.0 

0.0 

0.0 

7.8 

50.6 

42 

Q1928+7351 

0.982 

10.2 

0.0 

0.0 

0.0 

149.4 

1.9 

0.0 

0.0 

0.0 

12.4 

74.3 

45 

Q2243-1222 

0.966 

8.1 

0.3 

0.0 

0.1 

103.2 

0.0 

0.0 

2.2 

0.2 

8.1 

51.6 

49 

Q2344+0914 

0.963 

13.9 

0.8 

0.0 

0.0 

143.1 

1.8 

0.5 

1.6 

0.0 

5.3 

50.9 


Average 

0.962 

11.3 

0.9 

0.1 

0.0 

121.0 

0.9 

0.1 

1.0 

0.1 

8.4 

56.7 


STD 

0.017 

2.3 

0.9 

0.1 

0.0 

18.7 

0.8 

0.2 

0.9 

0.3 

1.9 

9.8 

Class II 

2 

Q0026+1259 

0.950 

9.6 

2.4 

0.1 

0.0 

92.1 

3.6 

0.1 

0.0 

0.0 

9.2 

19.4 

4 

Q0159-1147 

0.948 

2.8 

0.3 

0.9 

0.0 

52.8 

1.0 

0.0 

1.1 

0.8 

4.2 

16.8 

5 

Q0349-1438 

0.974 

1.8 

0.0 

0.0 

0.1 

54.9 

0.0 

0.0 

0.0 

0.0 

4.0 

28.0 

6 

Q0405-1219 

0.950 

6.9 

0.7 

1.1 

0.0 

92.5 

0.4 

0.0 

1.1 

0.0 

6.5 

35.0 

7 

Q0414-0601 

0.958 

6.1 

0.0 

0.1 

0.0 

114.0 

0.9 

0.0 

0.0 

0.0 

5.1 

40.1 

9 

Q0454-2203 

0.935 

11.6 

4.2 

0.0 

0.0 

105.7 

2.2 

0.3 

1.1 

0.0 

7.2 

33.9 

12 

Q0923+3915 

0.973 

7.2 

0.0 

0.1 

0.0 

95.1 

0.1 

0.0 

1.2 

0.0 

5.1 

46.4 

15 

Q0954+5537 

0.983 

3.3 

0.0 

0.4 

0.1 

42.6 

0.3 

1.1 

0.0 

0.0 

1.9 

17.9 

20 

Q1104+1644 

0.958 

5.7 

0.7 

0.0 

0.0 

104.5 

0.0 

0.0 

0.5 

0.0 

7.8 

54.1 

27 

Q1252+1157 

0.939 

6.1 

2.7 

2.4 

0.0 

72.7 

1.0 

0.2 

0.0 

0.5 

6.0 

25.5 

40 

Q1637+5726 

0.945 

4.5 

0.9 

0.2 

0.0 

77.6 

0.7 

0.0 

0.3 

0.0 

5.4 

33.2 

47 

Q2251+1552 

1.015 

7.1 

2.3 

0.0 

0.0 

60.1 

1.0 

0.0 

0.0 

0.0 

2.6 

23.2 


Average 

0.960 

6.1 

1.2 

0.4 

0.0 

80.4 

0.9 

0.1 

0.4 

0.1 

5.4 

31.1 


STD 

0.022 

2.8 

1.4 

0.7 

0.0 

23.7 

1.0 

0.3 

0.5 

0.2 

2.1 

11.6 

Class III 

16 

Q0959+6827 

0.865 

3.6 

4.0 

2.9 

0.1 

81.0 

3.1 

0.7 

1.3 

1.3 

12.1 

25.2 

17 

Q1001+2910 

0.949 

3.4 

3.4 

2.4 

0.0 

68.7 

2.0 

0.7 

2.9 

2.1 

9.3 

26.2 

21 

Q1115+4042 

0.957 

7.3 

4.3 

3.2 

0.0 

82.3 

1.8 

0.3 

2.5 

1.5 

8.6 

29.7 

23 

Q1148+5454 

0.928 

3.4 

3.9 

2.8 

0.0 

108.4 

1.1 

0.9 

1.8 

1.5 

14.1 

30.8 

26 

Q1248+4007 

0.937 

2.9 

3.0 

3.2 

0.1 

111.6 

2.1 

0.6 

1.6 

0.6 

12.2 

28.6 

28 

Q1259+5918 

0.954 

2.8 

2.5 

1.8 

0.0 

63.5 

1.0 

1.7 

2.7 

0.7 

10.1 

19.8 

29 

Q1317+2743 

0.963 

3.2 

0.8 

1.3 

0.0 

57.2 

1.5 

0.4 

0.5 

0.1 

7.1 

17.9 

30 

Q1320+2925 

0.977 

1.1 

3.9 

1.2 

0.0 

58.6 

1.0 

0.6 

0.4 

0.1 

5.7 

14.4 

36 

Q1444+4047 

0.925 

4.0 

3.0 

3.8 

0.1 

65.0 

3.6 

0.0 

0.3 

0.7 

9.9 

24.0 

43 

Q2145+0643 

0.943 

1.6 

0.0 

0.0 

0.0 

71.3 

0.6 

0.0 

0.0 

0.3 

1.6 

35.9 

44 

Q2201+3131 

0.928 

4.8 

0.6 

1.2 

0.0 

69.9 

0.0 

0.0 

0.0 

0.0 

7.6 

28.8 


Average 

0.939 

3.5 

2.7 

2.2 

0.0 

76.1 

1.6 

0.5 

1.3 

0.8 

8.9 

25.6 


STD 

0.029 

1.6 

1.5 

1.1 

0.0 

18.5 

1.1 

0.5 

1.1 

0.7 

3.5 

6.2 

Class IV 

3 

Q0044+0303 

0.965 

8.6 

1.3 

0.0 

0.0 

146.0 

1.5 

0.1 

1.7 

0.1 

11.3 

71.2 

10 

Q0624+6907 

0.888 

8.2 

5.1 

3.0 

0.0 

174.5 

3.9 

0.6 

2.5 

1.3 

11.3 

51.7 

19 

Q1100+7715 

0.948 

9.7 

1.7 

0.0 

0.0 

101.8 

1.0 

0.0 

0.0 

0.0 

8.5 

72.1 

31 

Q1322+6557 

0.925 

8.3 

1.8 

1.7 

0.0 

138.5 

2.7 

0.0 

3.3 

0.8 

9.3 

55.5 

33 

Q1402+2609 

0.865 

6.3 

3.1 

4.7 

0.0 

82.0 

1.8 

0.8 

2.3 

1.8 

5.1 

31.7 

35 

Q1427+4800 

0.944 

13.9 

1.8 

0.1 

0.0 

96.9 

0.7 

0.0 

0.6 

0.2 

16.9 

57.3 

38 

Q1544+4855 

0.860 

2.9 

5.7 

5.5 

0.0 

119.5 

1.0 

1.6 

5.1 

3.6 

16.5 

33.5 

46 

Q2251+1120 

0.937 

10.5 

0.8 

0.0 

0.1 

141.6 

4.5 

0.0 

2.1 

0.0 

7.8 

69.8 

50 

Q2352-3414 

0.937 

9.0 

0.8 

1.3 

0.0 

129.4 

0.4 

0.0 

0.7 

0.2 

6.6 

61.8 


Average 

0.919 

8.6 

2.4 

1.8 

0.0 

125.6 

1.9 

0.3 

2.0 

0.9 

10.4 

56.1 


STD 

0.038 

3.0 

1.8 

2.1 

0.0 

28.7 

1.5 

0.6 

1.6 

1.2 

4.1 

15.1 
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Table 2—Continued 


i 

Quasar 

< 5 F 

Ly / 3 /OVI 

Fe II 

Fe III 

C III * 

Lya 

N V 

Si II 

Si II 

C II 

Si IV 

C IV 




1025 

1071 

1123 

1176 

1216 

1240 

1263 

1307 

1335 

1397 

1549 




(A) 

(A) 

(A) 

(A) 

(A) 

(A) 

(A) 

(A) 

(A) 

(A) 

(A) 

Total 


Average 

0.947 

7.5 

1.6 

0.9 

0.0 

100.9 

1.2 

0.3 

1.0 

0.4 

8.1 

42.8 


STD 

0.031 

3.7 

1.5 

1.4 

0.0 

30.4 

1.1 

0.4 

1.1 

0.7 

3.3 

17.0 







